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DETECTION  PERFORMANCE  OF  OR-ING  DEVICE  WITH  PRE-  AND 
POST-AVERAGING:  PART  II  —  PHASE-INCOHERENT  SIGNAL 

INTRODUCTION 


BACKGROUND 

The  need  to  process  and  condense  large  amounts  of  data  is 
encountered  frequently  in  modern  Navy  systems  that  employ 
multiple  beams,  frequency  bins,  range  cells,  et  cetera.  One  way 
of  accomplishing  this  goal  is  by  or-ing  a  number  of  inputs  into  a 
single  output  —  that  is,  allowing  only  the  largest  of  a  set  of 
quantities  to  pass  on  for  further  processing  and  completely 
rejecting  the  remainder.  However,  since  this  or-ing  operation  is 
highly  nonlinear,  destroys  information,  and  tends  to  cause  small- 
signal  suppression,  some  pre-averaging  of  the  inputs  to  the 
or-ing  device  is  often  employed  in  an  effort  to  build  up  the 
signal-to-noise  ratio  (SNR)  prior  to  the  maximum  comparison. 
Furthermore,  at  the  or-ing  output,  some  additional  post-averaging 
is  frequently  used,  this  time  in  an  effort  to  build  up  the  SNR 
before  a  threshold  comparison  is  made  for  the  purpose  of 
declaring  a  signal  present  or  absent  (hypothesis  H^  versus  Hq, 
respectively) .  The  pertinent  block  diagram  is  displayed  in 
figure  1. 

There  are  N  channels  of  real  input  data  available  for 
processing,  namely,  (xn(t)}  for  1  <  n  <  N,  where  time  has  been 
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PRE¬ 
INPUTS  AVERAGER 


POST- 

OR-ING  AVERAGER  OUTPUT 


Figure  1.  Or-ing  with  Pre-  and  Post-Averaging 


normalized  so  that  time  sampling  instant  t  is  an  integer.  Under 
hypothesis  HQ,  there  is  Gaussian  noise  only  in  all  the  inputs, 
whereas  under  hypothesis  H^,  a  signal  is  also  present  in  one 
(unknown)  channel.  The  goal  of  the  processor  in  figure  1  is  to 
determine  signal  presence  with  a  high  detection  probability  P^, 
while  realizing  a  specified  acceptable  low  false  alarm 
probability  P^. 

Each  pre-averager  accumulates  K  statistically  independent 
consecutive  time  samples  of  its  corresponding  input  xn(t), 
yielding  output  yn(t),  which  is  then  subjected  to  or-ing  amongst 
N  competitors.  The  or-ing  output  at  time  t  is 

v(t)  -  maxfy^  t) ,  •  •  •  ,yN(  t) }  .  (1) 

Finally,  the  post-averager  accumulates  M  samples  of  its  input 
v(t)  and  compares  its  output  w(t)  with  a  fixed  threshold.  It  is 
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presumed  that  the  post-averager  waits  for  all  the  K  input  data 
samples  to  be  summed  and  or-ed  before  an  input  is  received;  then, 
another  block  of  K  input  samples  is  summed  and  or-ed,  leading  to 
the  next  input  to  the  post-averager.  That  is,  block  processing 
is  presumed. 


SCOPE 

This  technical  report  constitutes  part  II  of  a  series  of 
three  investigations  of  or-ing  with  pre-  and  post-averaging.  In 
particular,  three  different  input  signal  models  in  additive 
Gaussian  noise  are  of  interest: 

(I)  Random  (Gaussian)  signal, 

(II)  Phase-incoherent  signal,  and 

(III)  Coherent  (deterministic)  signal. 

Results  for  part  I,  the  Gaussian  signal,  have  been  documented  in 
reference  1;  the  reader  must  be  familiar  with  its  contents. 
Specifically,  the  details  on  block  and  overlap  processing,  as 
well  as  a  scaling  property  and  the  standard  of  comparison  for  the 
or-ing  system  in  figure  1,  are  presented  there,  along  with 
references  to  related  past  results. 
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The  current  report  addresses  part  II,  the  phase-incoherent 
signal  in  noise,  and  presents  numerous  receiver  operating 
characteristics  (ROCs)  that  completely  quantify  the  performance 
of  the  or-ing  system  in  figure  1.  For  all  three  signal  models 
above,  the  optimum  processors  have  been  determined  in  reference 
1,  appendix  A.  Thus,  the  exact  losses  that  the  or-ing  procedure 
causes  can  be  accurately  quantified. 

The  analytical  approach  utilized  here  can  be  extended  to 
include  quantizers  on  the  inputs  in  figure  1,  and  accurate  ROCs 
can  still  be  obtained  for  decision  variable  w(t).  In  fact,  there 
can  be  arbitrary  amounts  of  quantization  (L  levels),  and  the 
input  probability  density  functions  (PDFs)  can  be  completely 
arbitrary.  This  procedure  will  be  documented  in  part  III. 

For  block  processing  in  figure  1,  expressions  are  derived 
that  allow  for  accurate  evaluation  of  the  false  alarm  probability 
and  the  detection  probability  for  the  decision  variable  w(t). 
Furthermore,  this  is  accomplished  for  arbitrary  amounts  of 
pre-averaging  K,  arbitrary  amounts  of  or-ing  N,  arbitrary  amounts 
of  post-averaging  M,  and  arbitrary  input  SNRs.  No  approximations 
are  involved,  the  analysis  is  not  limited  to  mean  and  variance 
calculations,  and  no  appeal  is  made  to  the  central  limit  theorem. 
Rather,  the  approach  employs  a  judicious  combination  of  analysis 
with  computer-aided  numerical  calculations.  The  accuracy  of  the 
end  result  is  limited  only  by  the  accuracy  of  the  computer.  An 
entire  ROC  can  be  generated  in  minutes. 
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INPUT  DATA  MODELS 


In  this  section,  and  are  independent  zero-mean,  unit- 
variance  Gaussian  random  variables  (RVs),  and  e^  is  an  exponen¬ 
tial  RV  with  unit  mean,  that  is,  with  PDF  exp(-u)  for  u  >  0. 

PHASE-INCOHERENT  SIGNAL  IN  GAUSSIAN  NOISE 

Here,  in  part  II,  the  signal  channel  input  at  time  k  under  H^ 
is  an  envelope-squared  variate 

*k  -  l[(9k  +  Ak  cosek)2  +  (hk  ♦  Ak  slnek|2]  '  <2) 

where  signal  amplitude  A^  is  nonrandom  and  phase  0^  is  an  RV  with 
arbitrary  PDF,  that  is,  incoherent  phase.  The  input  mean  is 

*k  “  1  +  T  Ak  *  1  +  pk  '  (3) 

regardless  of  the  value  of  0^.  Since  x^  is  an  envelope-squared 

2 

quantity,  pk  =  A^/2  is  an  input  (per  sample)  power  SNR  measure. 

The  characteristic  function  (CF)  of  x^  under  H^ ,  conditioned 
on  a  particular  value  of  RV  0^,  is,  by  use  of  equation  (2), 

V'lV  •  exp(it%k)  -  j-rrt  exp(»k  r^n)  •  <4) 

But,  since  the  right-hand  side  of  CF  (4)  is  independent  of  0^, 
it  is  possible  to  use,  instead  of  form  (2),  the  data  value 
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xk  - §[%  +  Ak)2  +  hk] 


(5) 


for  the  signal  channel  input  at  time  k  under  H^.  On  the  other 
hand,  for  the  noise  channels,  Ak  is  zero  and  CF  (4)  reduces  to 
1/(1  -  i£),  which  corresponds  to  a  unit-mean  exponential  RV.  For 
the  noise  channel  inputs,  the  following  quantity  can  be  used: 


When  Ak  =  A  for  all  k,  that  is,  constant  signal  amplitudes 
with  time,  the  pre-averager  output  for  the  signal  channel  can  be 
expressed  as 


K 

y  -E 


k=l 


(g i  +  h?) 


K 

+  A  C 
k=l 


K  ^2- 
2 


=  y  +  A 


K 

e 

k=l 


K  - 
K  2 


(7) 


under  ,  where  y°  is  the  corresponding  noise-only  pre-averager 
output.  (However,  the  k-th  term  in  the  sum  for  y  cannot  be 
replaced  by  ek  +  A  gk  +  A  /2,  where  ek  is  an  independent 
exponential  RV,  because  the  CF  of  this  latter  combination  is 

1  A2 

1  ZTJI  expt  P  H(1  +  i^)]  ,  p  =  f  ,  (8) 


which  is  not  identical  in  form  to  equation  (4).  The  two  CFs 
do,  however,  have  the  same  mean  and  variance.)  The  ROCs  for 
pre-averager  (7)  can  be  parameterized  by  input  amplitude  SNR 
measure  A  or  by  input  power  SNR  measure  p. 
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Because  of  the  way  SNR  parameters  {p^}  appear  in  equation 
(4),  the  CF  of  the  pre-averager  output  can  be  found  in  the  closed 
form 


fy(S)  -  (1  il)  exP  (pK  i  -*ij  '  pK  ”  pk  '  (9) 

Therefore,  all  the  following  results  pertaining  to  the  phase- 
incoherent  signal  with  equal  SNRs  p  and  sums  of  squared  envelopes 
actually  hold  for  arbitrary  SNRs  { pk > ,  if  p  is  replaced  by  pR/K. 
However,  this  observation  is  of  limited  utility  unless  all  the  M 
sums  of  {pkJ,  one  for  each  post-averaging  interval,  are 
identical . 

OPTIMUM  PROCESSING  OF  AVAILABLE  DATA 

The  optimum  processor  for  signal  model  (2)  has  been  derived 
in  appendix  A  of  reference  1  for  the  case  where  channel  occupancy 
identification  is  required,  as  well  as  for  the  alternative  case 
where  channel  identification  is  irrelevant.  For  all  situations, 
the  optimum  processor  is  essentially  given  by  threshold 
comparison 

max  (c  x  (t)l  >  v  ;  (10) 

1  <n<N  lt=l  J  < 

see  equations  (A-13),  (A-19),  and  (A-26)  of  reference  1.  This 
operation  is  equivalent  to  performing  all  pre-averaging  and  no 
post-averaging  in  figure  1,  that  is,  taking  K  =  T,  the  total 
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time-bandwidth  product  available,  and  M  =  1.  Doing  so  defers  the 
nonlinearity  (or-ing)  until  after  all  possible  linear  combining 
(pre-averaging)  has  been  accomplished.  This  case  will  be 
thoroughly  evaluated  numerically  and  will  serve  as  the  basis  of 
comparison  for  the  various  combinations  of  pre— averaging  and 
post-averaging,  that  is,  K  <  T  or  M  >  1. 
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PERFORMANCE  ANALYSIS  FOR  BLOCK  PROCESSING 


Due  to  the  pre-averaging,  each  RV  yR(t)  in  figure  1  is  a  sum 
of  K  independent  RVs  with  identical  statistics.  Also,  all  N 
channel  inputs  are  statistically  independent  of  each  other. 

For  or-ing  to  be  present,  N  >  2  is  required  for  the  following 
calculations  and  considerations. 

For  a  noise-only  channel,  let  pQ  be  the  PDF  of  pre-averager 
output  Yn(t),  and  let  cQ  be  the  corresponding  cumulative 
distribution  function  (CDF),  namely,  cQ(u)  =  Prob{yn(t)  <  u|HQ}. 
For  a  signal  present  in  channel  1,  say,  let  p1  and  c1  be  the 
corresponding  PDF  and  CDF  of  y^t),  respectively. 

The  or-ing  output  v(t)  in  figure  1  is 

v ( t )  -  max{y1(t),...,yN(t)}  .  (11) 

Its  CDF  for  signal  present,  that  is,  hypothesis  ,  is 

cvl(u)  =  ci(u)  C (u ) ]N_1  ,  (12) 

which  leads  to  the  corresponding  PDF  of  v(t)  under  as 

pvl(u)  =  cvl(u)  =  pl(u)  ICqU)]"-1  +  c^u)  ( N— 1 )  [c0(u)]N_2  p0(u) 
=  (c0(u))N_2  [p1(u)  cQ(u)  +  ( N-l )  Pq ( u )  c-^u)]  .  (13) 

The  signal-absent  PDF  of  v(t)  under  HQ  follows  immediately  as 

pv0(u)  =  N  p0(u)  lc0(«)lN“1  •  (14) 
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The  corresponding  exceedance  distribution  functions  (EDFs), 
namely,  ey(u)  =  1  -  cv(u),  for  or-ing  output  v(t)  can  be 
expressed  in  terms  of  the  following  series  expansions,  which  are 
useful  for  very  small  EDF  values: 

1 

evl(u)  -  e1(u)  -  Cj^u)  )  )  [-e0(u)]n  ,  (15) 

ev0(u)  --E  C)  l-e0(u))“  .  (16) 

n=i 

The  CF  of  or-ing  output  v(t)  is  given  by  the  Fourier  transform 

fv(^)  =  exp(i£v(t) )  =  J  du  exp(i£u)  pv(u)  ,  (17) 

where  either  relevant  form,  pv^  or  pvq  above,  is  to  be  used. 
Finally,  for  block  processing  in  figure  1,  there  is  no  overlap  of 
the  data  used  in  the  post-averager,  meaning  that  the  processor 
output  w(t)  is  the  sum  of  M  independent  identically  distributed 
RVs .  Therefore,  the  CF  of  the  decision  variable  w(t)  in  figure  1 
is  simply 

fwU)  =  [fvU)lM  •  <18> 

One  final  Fourier  transform  is  required  to  manipulate  CF 
fw(  O  into  the  desired  EDF  of  processor  output  w(t);  see 
reference  2.  This  EDF  will  be  the  false  alarm  probability  Pf  or 
the  detection  probability  Pd  of  the  or-ing  processor,  depending 
on  whether  PDF  pvQ  or  pvl  is  used,  respectively. 


10 


The  fundamental  Fourier  transform  in  equation  (17)  will 
generally  have  to  be  accomplished  numerically  by  means  of  a  fast 
Fourier  transform  (FFT).  If  pre-averager  output  statistics  Pg 
and  Cg  can  be  found  in  closed  form  or  readily  computed  form,  a 
cascaded  FFT  approach  will  lead  to  exact  false  alarm  probability 
results  for  the  or-ing  processor  in  figure  1.  If  pre-averager 
output  probability  functions  p1  and  c1  can  also  be  readily 
evaluated,  accurate  detection  probability  results  can  be 
similarly  calculated  as  well.  Some  important  numerical 
considerations  for  this  cascaded  FFT  approach  are  presented  in 
appendix  B  of  reference  1. 

For  the  special  case  of  M  =  1  (no  post-averaging),  processor 
output  RV  w(t)  =  v(t)  in  figure  1,  and  the  corresponding  EDFs 
follow,  for  any  K,  N,  as 

ewl(u)  “  1  “  ci<u)  [c0(u)]N_1  ,  (19) 

ew0(u)  =  1  ~  tcQ(u)]N  .  (20) 

Now,  return  to  arbitrary  values  of  K,  N,  and  M,  and 
specialize  the  general  results  above  to  the  phase-incoherent 
signal  model  of  interest  here,  namely,  equations  (2)  and  (3). 
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PHASE-INCOHERENT  SIGNAL 


Define  p  -  K  p  and  use  (reference  3,  equation  6.5.11) 


,  ,  JSL,  xk 
e^(x>  =  i_.  rr  • 

k=0  K* 


n 


(21) 


For  constant  signal  amplitudes  and  u  >  0,  the  pertinent  PDFs  and 
CDFs  of  pre-averager  outputs  { yn ( t ) }  are  (using  equation  (9); 
reference  4,  equation  6.631  4;  and  reference  5,  equation  (1)) 

p0(u)  =  - — (Kfffi~u)  ,  cQ(u)  =  1  -  exp(-u)  eR_1(u)  ,  (22) 


K-l 


pl(u)  =  (p)  2  exP(-u"£)  ik_i(2(pu)J$)  ' 


(23) 


fx(5)  -  (1  -  U)  K  exp(p  ' 


(24) 


c1(u)  «  1  -  QR((2£)**,(2u)*)  ;  p  =  K  p  . 


(25) 


The  PDFs  for  or-ing  output  v(t),'  for  u  >  0,  then  follow  from 
equations  (14)  and  (13),  respectively,  as 

pv0(u)  =  N  - — (K-l'f l~Vl)  11  “  exP(-u)  eK-l(u)]N_1  '  (26) 

K-l 

pvi(u)  =  ( 1  -  exp(-u)  eK-1(u)  ]N_2  2  exp(-u-_p )  IR-1  (2 (£U ) **)  x 

K—l 

x  [1  -  exp(-u)  eR_1(u)]  +  (N-l)- — ,  J5  j  ~ui  ^1-Qr  (  ( 2p )  ** ,  ( 2u )  **  )  j  ]  . 

(27) 
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SPECIAL  CASES 


In  the  special  case  of  K  =  1,  that  is,  no  pre-averaging,  the 
CF  of  v(t)  with  or-ing  and  signal  present  is  available  in  closed 
form  by  Fourier  transforming  pv^(u)  to  obtain  (using  reference  4, 
equation  6.631  4,  and  reference  5,  equation  (24)) 


fvl<*> 


(-1) 


k+1 


(k  -  U)(k  +  1  -  U) 


exp(-p  iHhr^n)128’ 


for  any  N  > 


vl 


1  +  ♦  U  (V)  rfr)  •  <29> 


Finally,  when  p.  =  0,  equation  (28)  reduces  to  the  more  useful 
form  for  the  CF  as 

r  N  f  .  r.,-l  N  . 

fvo‘«  -  tni1  ■  *9] '  "vo  -  s  k  >  <3° 

this  alternative  form  was  derived  in  reference  1,  equation  (73). 


For  the  special  case  of  K  =  2  with  no  signal  present,  the  CF 
of  v(t)  is  available  by  Fourier  transforming  PDF  pvq(u),  yielding 


fv0<*>  " 


N! 


N-l 

E 

n=0 


(-1) 


n 


n 


(N-l-n) ! 


e 

k=0 


k+1 


(n-k)!  (n+l-i£) 


k+2 


(31) 


This  result  holds  for  arbitrary  amounts  of  or-ing,  i.e.,  N  >  1. 
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For  N  ■  2,  the  CF  of  v(t)  under  is  available  (after 
considerable  labor)  by  Fourier  transforming  PDF  Pvl(u)  to  obtain 


for  any  K  >  1,  where  -  1  -  i£  and  z2  =  2  -  i$.  For  K  -  2,  the 
mean  follows  as  "vl  *  2(1  +  p)  +  exp(-p)  (6  +  p)/8,  using  £  =  2p. 


This  form  for  CF  fvQ(S)  in  equation  (33)  is  not  obtained  by 
setting  £  =  0  in  equation  (32);  however,  numerical  investigation 
of  both  forms  reveals  that  they  are  equal.  For  K  -  2,  the  CF  in 
equation  (33)  further  specializes  to 


fv0<*>  " 


(35) 
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The  CF  of  decision  variable  w(t)  for  signal  present  is  given 
by 

fKlU)  -  lfvlU)]M  ■  (36) 

When  CF  fw^(£)  is  available,  a  Fourier  transform  will  yield  the 
detection  (or  false  alarm)  probabilities  for  any  K,  N,  and  M. 

When  M  =  1,  the  EDFs  in  equations  (19)  and  (20)  can  be  combined 
with  results  (22)  and  (25)  to  immediately  yield  closed  form 
expressions  for  P^  and  Pf  for  any  K,  N,  and  SNR  p. 

When  N  ■  1  (no  or-ing)  and  K  and  M  are  arbitrary,  the  results 
simplify  to  output  PDF  and  EDF 

T-l 

pwl(u)  "  (b)  2  exP(“u“B)  jt-1  (2(Bu),S)  for  u  >  0  '  (37) 

00 

Pd(v)  "  J  du  pwl(u)  =  QT((2B)Js,  (2V)55]  for  v  >  0  ,  (38) 

v 

respectively,  where  T  -  KM  and  B  -  Tp  -  KMp.  Obviously,  for 
N  -  1  here,  only  the  product  KM  matters;  this  presumes  all  unity 
weights  in  both  the  pre-  and  post-averagers. 


NUMERICAL  APPROACH 

When  the  CF  fv(£)  of  or-ing  output  v(t)  can  be  derived  in 
closed  form,  such  as  in  special  cases  (28)  through  (35),  f  ($) 
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can  be  used  directly  in  equation  (18)  to  find  the  CF  of  decision 
variable  w(t).  This  route  is  preferred  because  CF  fv(£)  can 
decay  rather  slowly  with  $  for  some  values  of  K  and  N.  In 
particular,  if  the  properties 

1  -  exp(-u)  eR_^(u)  «  uK  as  u  ->  0+  ,  (39) 

In(x)  «  xn  as  x  •+  0+  ,  (40) 

1  -  Q_(a,x)  «  x2K  as  x  -*  0+  (41) 

x\ 

are  used,  it  follows  that  both  PDFs  pvg(u)  and  pvl(u)  in 

equations  (26)  and  (27)  are  proportional  to  uKN-^  as  u  -»  0+.  For 

small  KN,  this  abrupt  change  (from  the  constant  value  0  for 

KN 

u  <  0)  leads  to  a  slow  decay  of  fv(0  with  $,  namely,  1/S  ,  with 

the  attendant  aliasing  problems;  see  reference  1.  Therefore, 
closed  form  results  for  CF  fv(S),  when  available,  can  be  very 
useful  in  terms  of  avoiding  some  numerical  problems  in  the 
discrete  Fourier  transform  of  pv(u).  When  this  is  not  possible, 
the  closed-form  results  for  the  or-ing  output  PDFs  given  in 
equations  (26)  and  (27)  are  used,  instead,  in  the  cascaded  FFT 
procedure  outlined  in  equations  (17)  and  (18). 

The  very  special  case  of  K  ■  1,  B  ■  1  can  be  solved 
explicitly  for  P^  and  P^,  as  shown  in  appendix  A.  Also,  the 
required  input  SNRs  for  two  typical  operating  points  are 
tabulated  in  this  appendix.  For  other  K , M  values,  the 
alias-suppression  technique  of  appendix  B  has  been  used  to 
evaluate  the  ROCs  directly  from  PDF  equations  (26)  and  (27). 
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QUANTITATIVE  PERFORMANCE  RESULTS 

This  section  gives  specific  quantitative  detectability 
results  for  or-ing  with  various  amounts  of  pre-averaging  K  and 
post-averaging  M.  In  particular,  numerous  ROCs  are  presented  in 
appendixes  C,  D,  E,  and  F  for  the  cases  of  KM  =  4,  16,  64,  and 
256,  respectively.  The  particular  cases  run  are  indicated  in 
table  1;  a  total  of  124  ROCs  are  represented.  For  N  =  1,  since 
only  the  product  KM  matters,  only  one  ROC  need  be  presented  (the 
first  curve  in  each  appendix)  under  the  title  labeled  M  =  1. 


Table  1.  Cases  Run 


K 

M 

N 

4 

1  • 

1,2,4,8,16,32 

Appendix  C 

2 

2 

1+15  ROCs 

1 

4 

16 

1 

1,2,4,8,16,32 

Appendix  D 

8 

2 

1+25  ROCs 

4 

4 

2 

8 

1 

16 

64 

1 

1,2,4,8,16,32 

Appendix  E 

32 

2 

1+35  ROCs 

16 

4 

8 

8 

4 

16 

2 

32 

1 

64 

256 

1 

1,2,4,8,16,32 

Appendix  F 

128 

2 

1+45  ROCs 

64 

4 

32 

8 

16 

16 

8 

32 

4 

64 

2 

128 

1 

256 
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REQUIRED  INPUT  SIGNAL-TO-NOISE  RATIOS 


The  numerous  ROCs  in  appendixes  C  through  F  enable  a  user  to 
easily  investigate  and  determine  the  amount  of  system  input  SNR 
required  for  the  or-ing  device  in  figure  1  to  achieve  a  specified 
level  of  detectability  performance,  in  terms  of  a  desired  false 
alarm  probability  and  detection  probability  P^,  over  a  wide 
range  of  parameter  values.  To  partially  condense  this  voluminous 
information,  a  standard  operating  point  (SOP)  is  defined  here  as 
P^  =  IE-3  and  P^  =  0.5,  and  a  high-quality  operating  point  (HOP) 
is  defined  as  Pf  =  IE-6  and  Pd  =  0.9.  The  corresponding  required 
values  of  the  input  SNR  p(dB)  =  10  log(p),  as  determined  from  the 
ROCs  in  appendixes  C  through  F,  are  listed  in  tables  2  through  9. 
(These  accurate  values  were  interpolated  from  the  ROCs  while  the 
false  alarm  and  detection  probability  numbers,  P^  and  P^,  were 
still  resident  in  the  computer;  eyeball  interpolation  from  the 
plotted  ROCs  cannot  be  accomplished  this  accurately. ) 

The  input  SNR  requirements  are  plotted  in  figures  2  through 
9.  As  expected,  the  required  input  SNR  p  increases  monotonically 
with  the  number  (N)  of  channels  or-ed  if  K  and  M  are  held  fixed. 
Also,  the  required  input  SNR  p  increases  monotonically  with  the 
amount  (M)  of  post-averaging  employed  if  N  and  product  KM  are 
held  fixed;  alternatively,  the  required  input  SNR  p  decreases 
with  the  amount  (K)  of  pre-averaging  employed  if  N  and  product  KM 
are  held  fixed.  The  exact  rates  can  be  determined  from  figures  2 
through  9  and  the  ROCs  in  the  appendixes. 
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Table  2.  Required  Input  SNR  p(dB)  for  =  IE-3,  P.  =  0.5, 
KM  =  4,  Phase-Incoherent  Signal 


K 

M 

| 

m 

ED 

N  *  8 

N  =  16 

N  =  32 

4 

1 

B 

4.15 

4.50 

4.82 

5.10 

5.38 

2 

2 

4.40 

4.95 

5.43 

5.87 

6.27 

1 

4 

4.62 

5.40 

6.07 

6.68 

7.21 

Table  3.  Required  Input  SNR  p(dB)  for  P^  =  IE-6,  P.  =  0.9, 
KM  s  4,  Phase-Incoherent  Signal 


K 

i 

M 

N  -  1 

N  =  2 

N  =  4 

N  -  8 

N  -  16 

N  =  32 

4 

1 

8.25 

8.40 

8.54 

8.68 

8.83 

8.95 

2 

2 

8.25 

8.53 

8.81 

9.05 

9.29 

9.50 

1 

4 

8.25 

8.71 

9.14 

9.53 

9.90 

10.23 

Table  4.  Required  Input  SNR  p(dB)  for  Pf  =  IE-3,  P,  =  0.5, 
KM  =  16,  Phase-Incoherent  Signal 


K 

M 

N  =  1 

N  =  2 

N  =  4 

00 

It 

a 

N  =  16 

N  =  32 

16 

i 

ma 

0.25 

0.55 

0.83 

1.08 

1.32 

8 

2 

wBSSm 

0.48 

0.97 

1.41 

1.79 

2.12 

4 

4 

-0.08 

0.74 

1.44 

2.03 

2.56 

3.02 

2 

8 

-0.08 

0.96 

1.88 

2.68 

3.37 

3.97 

1 

16 

-0.08 

1.11 

2.25 

3.26 

4.14 

4.92 

Table  5.  Required  Input  SNR  p(dB)  for  Pf  -  IE-6 ,  P.  =  0.9, 
KM  =  16,  Phase-Incoherent  Signal 


K 

M 

N  =  1 

N  =  2 

■n 

00 

It 

£ 

N  =  16 

N  -  32 

16 

1 

3.84 

3.97 

4.09 

4.22 

4.34 

4.45 

8 

2 

3.84 

4.09 

4.33 

4.54 

4.75 

4.94 

4 

4 

3.84 

4.28 

4.66 

5.02 

5.34 

5.63 

2 

8 

3.84 

4.48 

5.07 

5.59 

6.04 

6.45 

1 

16 

3.84 

4.69 

5.48 

6.18 

6.80 

7.36 
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Table  6.  Required  Input  SNR  p(dB)  for  P^  =  IE-3,  =  0.5, 

KM  =  64,  Phase-Incoherent  Signal 


K 

M 

ID 

N  =  4 

N  =  8 

VO 

1—1 

II 

£ 

N  -  32 

64 

1 

-3.59 

-3.28 

-3.01 

-2.75 

-2.53 

-2.32 

32 

2 

-3.59 

-3.06 

-2.61 

-2.22 

-1.88 

-1.57 

16 

4 

-3.59 

-2.82 

-2.16 

-1.61 

-1.14 

-0.73 

8 

8 

-3.59 

-2.57 

-1.71 

-1.00 

-0.39 

0.14 

4 

16 

-3.59 

-2.37 

-1.30 

-0.40 

0.39 

1.05 

2 

32 

-3.59 

-2.22 

-0.96 

0.15 

1.12 

1.95 

1 

64 

-3.59 

-2.14 

-0.72 

0.58 

1.77 

2.79 

Table  7.  Required  Input  SNR  p(dB)  for  P^  =  IE-6,  P^  =  0.9, 
KM  =  64,  Phase-Incoherent  Signal 


K 

M 

ESI 

N  =  4 

N  =  8 

N  -  16 

N  =  32 

64 

i 

-0.08 

0.03 

0.15 

0.26 

0.37 

0.46 

32 

2 

-0.08 

0.14 

0.36 

0.55 

0.74 

0.91 

16 

4 

-0.08 

0.33 

0.68 

0.99 

1.28 

1.53 

8 

8 

-0.08 

0.54 

1.07 

1.53 

1.94 

2.29 

4 

16 

-0.08 

0.77 

1.50 

2.13 

2.67 

3.15 

2 

32 

-0.08 

0.97 

1.92 

2.74 

3.44 

4.05 

1 

64 

-0.08 

1.11 

2.26 

3.28 

4.17 

4.96 

Table  8.  Required  Input  SNR  p(dB)  for  P^  =  IE-3,  P.  =  0.5, 
KM  =  256,  Phase-IncoherentrSignal 


K 

M 

N  =  1 

N  -  2 

N  =  4 

N  =  8 

N  =  16 

N  =  32 

256 

i 

-6.87 

-6.57 

-6.31 

-6.08 

-5.86 

-5.65 

128 

2 

-6.87 

-6.36 

-5.94 

-5.56 

-5.24 

-4.96 

64 

4 

-6.87 

-6.11 

-5.50 

-4.99 

-4.54 

-4.16 

32 

8 

-6.87 

-5.87 

-5.06 

-4.38 

-3.81 

-3.33 

16 

16 

-6.87 

-5.65 

-4.63 

-3.78 

-3.07 

-2.47 

8 

32 

-6.87 

-5.48 

-4.25 

-3.21 

-2.35 

-1.61 

4 

64 

-6.87 

-5.36 

-3.96 

-2.72 

-1.67 

-0.78 

2  ' 

128 

-6.87 

-5.29 

-3.75 

-2.34 

-1.09 

-0.02 

1 

256 

-6.87 

-5.28 

-3.67 

-2.10 

-0.64 

0.64 
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Table  9.  Required  Input  SNR  p(dB)  for  P^  =  IE-6,  =  0.9, 

km  ■  256,  Phase-Incoherent  signal 


K 

M 

N  -  1 

N  =  2 

2 

II 

00 

N  =  16 

N  -  32 

256 

1 

-3.62 

-3.51 

-3.41 

-3.31 

-3.21 

-3.12 

128 

2 

-3.62 

-3.41 

-3.21 

-3.04 

-2.87 

-2.71 

64 

4 

-3.62 

-3.25 

-2.92 

-2.62 

-2.37 

-2.13 

32 

8 

-3.62 

-3.03 

-2.54 

-2.11 

-1.75 

-1.43 

16 

16 

-3.62 

-2.81 

-2.13 

-1.54 

-1.05 

-0.63 

8 

32 

-3.62 

-2.58 

-1.69 

-0.96 

-0.33 

0.21 

4 

64 

-3.62 

-2.40 

-1.31 

-0.38 

0.42 

1.09 

2 

128 

-3.62 

-2.26 

-0.99 

0.14 

1.11 

1.96 

1 

256 

-3.62 

-2.19 

-0.78 

0.54 

1.73 

2.76 

The  actual  numerical  procedure  used  to  generate  the  ROCs  in 
appendixes  C  through  F  is  based  upon  appendix  B  and,  in 
particular,  on  smoother  function  b^(u)  and  equation  (B-14). 
However,  an  alternative  approach,  presented  in  appendix  G,  for 
evaluation  of  a  CF  directly  from  samples  of  a  CDF  utilizes  a 
model  CDF  function  that  mimics  the  general  behavior  of  the  CDF  of 
interest.  Proper  control  of  aliasing  and  efficient  use  of  an  FFT 
are  demonstrated  in  appendix  G,  along  with  an  example  of  the 
application  of  the  alternative  procedure  to  a  Gaussian  CF.  The 
inadequacy  of  simply  using  CDF  differences  as  estimates  of  the 
PDF  is  also  illustrated  by  example  there. 
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TRADEOFF  BETWEEN  PRE-  AND  POST-AVERAGING 


For  a  fixed  amount  of  total  time— bandwidth  product,  that  is, 
KM  constant,  it  is  possible  to  shift  some  of  the  pre-averaging  to 
post-averaging,  or  vice-versa.  For  example,  if  KM  *  16,  this 
could  be  accomplished  by  taking  K  =  16,  M  =  1,  which  corresponds 
to  no  post-averaging,  and  switching  to  K  =  1,  M  =  16,  which 
corresponds  to  no  pre-averaging.  Alternatively,  intermediate 
values  such  as  K  =  4,  M  =  4,  could  be  used,  still  keeping 
KM  =  16. 


In  all  cases,  increasing  M  at  the  sake  of  K  requires  larger 
input  SNRs  to  maintain  the  same  operating  point;  see  figures  2 
through  9.  This  is  related  to  the  fact  that  the  optimum 
processor  conducts  all  pre-averaging  before  any  nonlinear 
operations;  see  appendix  A.  For  example,  to  maintain  the  SOP 
when  KM  =  16,  table  4  indicates  that  increases  of  0.86  dB,  1.70 
dB,  2.43  dB,  3.06  dB,  and  3.60  dB  are  required  for  N  =  2,  4,  8, 
16,  and  32,  respectively,  if  the  switch  is  made  from  all  pre¬ 
averaging  to  all  post-averaging.  Alternatively,  for  the  HOP,  the 
corresponding  differences  in  table  5  are  smaller  (0.72  dB,  1.39 
dB,  1.96  dB,  2.46  dB,  and  2.91  dB  for  N  =  2,  4,  8,  16,  and  32, 
respectively) . 

For  a  HOP  such  as  Pf  =  IE-6  and  Pd  =  0.9,  larger  input 
SNRs  are  naturally  required  to  achieve  the  higher  level  of 
performance.  At  these  larger  input  SNRs,  the  or-ing  process 
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is  more  frequently  dominated  by  the  signal-bearing  channel. 
Therefore,  as  the  number  (N)  of  channels  or-ed  increases,  the 
required  increases  in  SNR  are  less  severe  than  for  a  lower- 
quality  operating  point,  such  as  Pf  =  lE-3  and  Pd  =  0.5.  For 
example,  for  K  =  16  and  M  =  1,  increasing  N  from  2  to  32  requires 
a  1.07-dB  SNR  increase  to  maintain  the  SOP,  whereas  only  an 
additional  0.48  dB  is  needed  at  the  HOP. 

A  related  observation  is  that,  for  a  given  configuration 
(fixed  K  and  M),  the  increase  in  SNR  required  to  maintain  the 
SOP,  as  the  amount  of  or-ing  increases,  is  more  severe  for  the 
larger  M  values.  Thus,  for  KM  -  16,  as  N  increases  from  2  to  32, 
a  1.07-dB  SNR  increase  (1.32  -  0.25)  suffices  to  maintain  the  SOP 
for  M  =  1,  whereas  a  3.81-dB  increase  (4.92  -  1.11)  is  required 
for  H  ■  16.  The  corresponding  increases  at  the  HOP  are  0.48  dB 
(4.45  -  3.97)  for  M  *  1  versus  2.67  dB  (7.36  -  4.69)  for  M  =  16, 
as  N  increases  from  2  to  32.  The  behaviors  are  similar  for  the 
other  values  of  KM,  although  the  required  input  SNR  p  values  are 
smaller  as  KM  increases,  due  to  the  additional  observation  times 
allowed. 
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Figure  9.  Required  Input  SNR  for  P£  =  IE- 6,  Pd  =  0.9,  KM  =  256 

Phase-Incoherent  Signal 
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SUMMARY 


For  a  phase-incoherent  signal  in  Gaussian  noise,  closed  forms 
have  been  derived  for  the  PDFs  at  the  output  of  a  pre-averager 
followed  by  an  or-ing  operation.  These  forms  have  been 
numerically  doubly  Fourier-transformed  to  yield  accurate  results 
for  the  EDFs  at  the  output  of  a  post-averager,  under  both 
hypotheses  Hq  and  H^,  for  numerous  values  of  the  parameters  of 
the  complete  processor  (figure  1).  Such  results  enable 
investigation  of  the  or-ing  processor  for  false  alarm 
probabilities  in  the  range  of  IE-6  and  smaller;  there  is  no 
need  to  resort  to  lengthy  simulations. 

Numerous  ROCs  have  been  generated,  which  enable  a  user  or 
system  designer  to  quickly  assess  the  losses  to  be  expected  from 
employing  or-ing  in  the  processor.  Also,  quantitative  evaluation 
of  the  tradeoffs  between  pre-averaging  and  post-averaging  has 
been  conducted.  Finally,  since  the  enclosed  tabulations  will 
undoubtedly  not  cover  all  cases  of  practical  interest,  a  MATLAB 
program  for  evaluation  of  additional  ROCs  is  presented  in 
appendix  H. 

The  last  case  of  a  deterministic  signal  in  additive  Gaussian 
noise  is  currently  under  investigation.  Tabulation  of  the 
ROCs  for  part  III  and  tradeoffs  between  pre-  and  post-averaging 
will  be  conducted  in  a  similar  manner  and  for  the  same  parameter 
values  as  accomplished  here  and  in  reference  1. 
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APPENDIX  A  -  REQUIRED  INPUT  SNR  FOR  K  =  1,  M  =  1 

The  cases  of  KM  -  4,  16,  64,  and  256  for  the  phase-incoherent 
signal  have  been  investigated  elsewhere  in  this  report.  In  this 
appendix,  to  complete  the  sequence,  results  for  KM  =  1  are 
presented.  Also,  complete  analytic  results  are  attainable  in 
this  very  special  case,  allowing  for  direct  calculation  of  the 
required  input  SNRs  to  realize  a  specified  level  of  performance 
(Pf  and  Pd)  at  the  output  of  the  or-ing  processor  in  figure  1. 

Use  of  equations  (19)  and  (20)  with  M  «  1  and  any  K  results 


in 

Pd  =  1  -  Cl(u)  [ c0(u)  ]N_;L  ,  ( A— 1 ) 

Pf  “  1  -  £c0(u)jN  •  (A- 2 ) 

Then,  in  particular,  for  the  phase-incoherent  signal  and  K  =  1, 
equations  (22)  and  (25)  yield,  for  any  N  and  u  >  0, 

c0(u)  =  1  -  exp(-u)  ,  C;L(u)  =  1  -  q((2p)Js,(2u)J5)  ,  (A-3) 

Pd  =  1  -  [l  -  Q((2p)Js,(2u)Js]]  [1  -  exp(-u)  ]N_1  ,  (A-4) 

Pf  =  1  -  [ 1  -  exp(-u)]N  .  (A-5) 

The  ROCs  for  the  phase-incoherent  signal  for  K  ■  1,  M  =  1,  and 


any  N  and  p  could  be  generated  directly  from  these  last  two 
relations,  if  desired. 
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Define  the  auxiliary  parameters 


a 


(1  - 


Pf> 


1/N 


( A-6 ) 


which  can  be  computed  once  the  desired  and  are  specified. 
Then,  the  threshold  u  required  to  realize  a  specified  P^  follows 
from  equations  (A-5)  and  (A-6)  as 


u  =  — In ( 1  —  a)  . 


(A— 7 ) 


In  a  similar  vein,  the  required  input  SNR  p  follows  from  equation 
(A-4),  using  the  threshold  u  given  by  equation  (A-7),  as  the 
solution  of  the  transcendental  equation 

<z[(2p)hA-2  ln(  1  -  a))*5]  =  1  -  0  a  .  (A-8) 

Results  for  the  SOP  and  HOP  are  shown  in  tables  A-l  and  A-2 . 


Table  A-l.  Required  Input  SNR  p(dB)  for  P^  =  IE-3,  P^  =  0.5, 
K  =  1,  M  =  1,  Phase-Incoherent  Signal 


N 

1 

2 

4 

8 

16 

32 

p(dB) 

8.06 

8.51 

8.91 

9.28 

9.62 

9.94 

Table  A-2.  Required  Input  SNR  p(dB)  for  Pf  ■  IE-6, 
K  =  1,  M  =  1,  Phase-Incoherent  Signal 

Fd  -  0.9, 

N  1  2  4  8  16 

32 

P(dB)  13.18  13.36  13.53  13.69  13.84 

13.99 
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APPENDIX  B  -  SUPPRESSION  OF  ALIASING  EFFECTS  NEAR  ORIGIN 


In  this  appendix,  it  is  assumed  that  the  PDF  p(u)  and  the  CDF 
c(u)  of  an  RV  are  given  and  that  it  is  desired  to  evaluate  an 
accurate  version  of  the  corresponding  CF  f(5),  especially  near 
the  origin  5=0.  This  is  very  important  when  relations  (18)  and 
(36)  are  used  for  large  M,  where  the  near-origin  behavior  of  CF 
fv(5)  is  the  most  critical  feature. 


To  this  end,  define  a  unit-area  localized  function  b(u), 
centered  at  u  =  0.  Several  examples  are  displayed  in  figure  B-l. 


bQ(u) 

bx(u) 

b2(u) 

UNIT-  J 

1/A  / 

V 1/A 

AREA 

/ 

\ 

IMPULSE 

/ 

\ 

u 

u  / 

\  u 

0  A 

2 

l 

d 

\ 

A 

Figure  B-l.  Possible  Localized  Functions 


Also,  define  convolution 


s(u)  =  p(u)  ©  b(u)  =  J  dt  p(t)  b(u-t)  =  J  dt  c(t)  b'(u-t)  .  (B-l) 


Then,  for  the  bg(u)  function,  there  follows  convolution 
Sq(u)  =  p(u);  for  the  b^(u)  function,  the  convolution  is 

u+A/2 

sx(u)  -  ^  J  dt  p(t)  =  i[c(u  +  j]  -  c  (u  -  jj]  ;  ( B-2 ) 

u-A/2 


B-l 


and  for  the  b2<u)  function,  the  convolution  is 

u+A  A 

s2(u)  -  i  j  dt  p(t)  (l  -  ,  _1 


u-A 


dt  [c(u+t)  -  c(u-t) ] . (B-3) 


Finally,  define  an  impulsive  approximation  to  PDF  p(u)  as 

p(u)  S  S  an  8 (u  -  nA)  ,  (B-4) 

n 

where  impulse  area  an,  defined  as 

an  e  a  s(na)  =  a  J  dt  p(t)  b(na  -  t)  ,  (B-5) 

is  a  localized  average  of  the  PDF  p(u),  centered  at  nA.  For  the 
b^(u)  example,  the  area  is 

aln  =  A  si<nA)  =  c[(n  +  Jj)A]  -  c[(n  -  *s)A]  ,  ( B— 6 ) 

and  for  b^{u) ,  the  area  is  aQn  =  A  Sg(nA)  =  A  p(nA). 

The  impulsive  PDF  in  equation  (B-4)  is  expressible  as 

p(u)  =  £  A  s(nA)  8(u  -  nA)  =  A  s(u)  £  8(u  -  nA)  s  s(u)  AS.(u)  , 
n  n 

( B-7 ) 

where  8^(u)  is  an  infinite  impulse  train  of  spacing  A  in  u.  The 
corresponding  CF  is 

?(£)  s  T  du  exp(iSu)  p(u)  =  £  a  exp(i^nA)  - 
J  n 

.  r  au  exp(i£u)  s(u)  A8^(u)  - 
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( B-8 ) 


-  3?  ■(«)  •  [2"  s2B/i<£>]  -  S  «(«-  «H)  > 

where  function 

SU)  s  J  du  exp(iSu)  s(u)  =  J  du  exp(i$u)  [p(u)  ©  b(u)]  = 

«  f<S)  B(S)  ,  ( b-9 ) 

with  definition 

B(£)  ■  J  du  exp(iSu)  b(u)  .  (B-10) 

Therefore,  the  impulsive-approximation  CF  is  given  by 

*<«  -  S  *(«  -  •(«  -  - 

■  f(S)  B<£)  +  E  f(*  -  b[?  -  n^|j  .  ( B-ll ) 


The  three  earlier  examples  of  b(u)  in  figure  B-l  can  be 
Fourier  transformed  according  to  equation  (B-10),  yielding 


B0U)  =  1  , 


B1(5) 


sin( A£/2 ) 
A 1/2 


B,u>  .  fsiniMZlil2 
2^’  L  A£/2  J 


( B-12 ) 


Thus,  for  no  smoothing  bQ(u),  the  standard  result  is 


f0u> 


S  a0n  exP( i )  =  4  JJ  p(nA)  exp(HnA)  = 
n  n 


-  fU>  +  2  f(*  -  n2i)  , 

n^n  v  a' 


( B-l 3 ) 


where  the  last  term  represents  the  sum  of  all  the  undesired 
aliasing  lobes  in  the  £  domain,  if  pdf  p(u)  is  sufficiently 
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smooth  for  all  u,  such  that  CF  f(5)  decays  very  rapidly  with  5, 
then  sampling  increment  A  in  equation  ( B— 4 )  can  easily  be  chosen 
small  enough  that  an  essentially  alias-free  region  near  the 
origin  5  =  0  can  be  found  where  CF  approximation  2Q(5)  is 
essentially  equal  to  the  desired  CF  f(5). 

However ,  if  PDF  p(u)  is  not  sufficiently  smooth  (for  example, 
u  exp(-u)  for  u  >  0),  then  CF  f(5)  decays  slowly,  causing 
significant  aliasing  effects,  due  to  the  last  term  in  equation 
( b-13 )  ♦  These  effects  can  be  suppressed  by  using  first-order 
smoothing  b^(u),  giving  an  alternative  approximate  CF  that 
utilizes  CDF  samples,  namely, 

2- (5)  =  Y,  ain  exp(i£nA)  =  £  [c(nA+J$A) 

1  n  n 

=  f<o  b  m  +  Y  f (s  -  Bi(s 

But,  from  equation  ( B— 12 ) ,  since 

B^^n/A)  =  sin(nrc)/(nn)  =  0  for  n  ^  0  ,  (B-15) 

all  the  aliasing  lobe  contributions  (real  and  imaginary)  are  zero 
at  5  -  0,  the  center  of  the  desired  lobe  of  equation  (B-14).  In 
fact,  the  real  part  of  the  aliasing  components  of  22(5)  also  has 
zero  slope  at  5  =  0;  however,  the  imaginary  part  has  a  nonzero 
slope  at  5  -  0.  These  effects  result  in  a  small  aliasing 
contribution  in  the  immediate  neighborhood  of  5  =  0,  leading  to 
the  approximation 


-  c(nA-SjA)]  exp(i£nA)  = 

-  n^£)  .  (B-14) 


B— 4 


=  fU)  B1(^)  for  l  near  0  . 


( B-16 ) 


Finally,  since  B^(£)  is  a  known  function,  the  desired  CF  f(£) 
can  be  obtained  approximately  as 


fU) 


^U) 

b^TT) 


for  £  near  0  . 


( B-17 ) 


If  CF  f(£)  must  be  raised  to  a  power,  as  in  the  averaging  of 
several  RVs,  each  with  PDF  p(u),  this  approach  is  particularly 
attractive  for  circumventing  aliasing  effects;  the  larger  the 
power,  the  more  useful  is  this  alternative.  Notice  that  the 
division  by  B^(£)  in  equation  (B-17)  must  be  done  before  raising 
the  CF  to  the  power  of  interest. 


For  the  smoother  bjfu)  in  figure  B-l,  there  follows,  from 
equation  (B-12), 


B2(n2n/A)  =  0  and  B£(n2jt/A)  =  0  for  n  ^  0  .  (B-18) 

Thus,  the  total  aliasing  component  of  i^tS)  and  *ts  sl°Pe 
(real  and  imaginary)  are  zero  at  £  =  0.  This  leads  to  the 
approximation  for  the  desired  CF  as 


fU) 


-  J2(S) 
*  B2(t) 


for  £  near  0  . 


( B-19 ) 


The  utility  of  equation  (B-19)  depends  on  the  ability  to  evaluate 
&2(nA)  in  equation  (B-3)  simply. 
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Equations  (B-6)  and  (B-14)  presume  that  samples  of  the  CDF 
c ( u )  are  available  at  in-between  points  u  -  (n+*s)A.  If  the  CDF 
samples  are  only  available  at  the  points  u  =  nA,  a  modified 
smoothing  function  b3(u)  can  be  employed.  The  pertinent 
equations  are 


b3(u) 


pL/(  2A)  for  |u|  <  A'j 
L  0  otherwise  J 


B3(5)  =  Sin^^}-  ,  ( B-20 ) 


u+A 


s3(u)  =  dt  p(t)  ^  = 


2Atc(u+A) 


-  C(u-A) ]  , 


( B-21 ) 


u-A 


a3n  -  A  s3(nA)  =  ^-[c(nA+A)  -  c(nA-A)]  , 


( B-22 ) 


?,(£)  -  £  a3n  exp(i^nA)  =  [c(nA+A)  -  c(nA-A)]  exp(i^nA)  - 


n 


n 


fm  B  U)  +  £  f  U  -  b3(?  -  n^)  , 

n^O  v  '  v 


( B-23 ) 


2qU) 

f(U  s  bJuT  for  *  near  0  * 


( B— 24 ) 


Since  B3U)  in  equation  (B-20)  is  zero  at  |S|  «  n/A,  it  is 
necessary  to  keep  |£|  <  n/A  when  equation  (B-24)  is  used;  this 
limit  is  twice  as  restrictive  as  the  one  for  the  use  of  B^(£)  in 
equations  (B-12)  and  (B-17).  Relations  (B-23)  and  (B-24)  afford 
an  alternative  method  of  proceeding  directly  from  samples  of  the 
CDF  c(u)  to  the  CF  f(£).  The  quantity  on  the  right-hand  side  of 
equation  (B-21)  can  be  interpreted  as  an  estimate  of  the  slope  of 
the  CDF  c(u)  at  u,  that  is,  an  estimate  of  the  PDF  p(u)  at  u. 
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APPENDIX  C  -  ROCS  FOR  KM  =  4,  PHASE- INCOHERENT  SIGNAL 


This  appendix  contains  the  ROCs  for  or-ing  with  pre-  and 
post-averaging  when  the  time-bandwidth  product  KM  is  fixed  at  4; 
the  possible  combinations  (from  table  1)  are  repeated  here: 


K 

M 

N 

4 

1 

1,2,4,8,16,32 

2 

2 

1 

4 

For  N  =  1,  only  the  product  KM  matters;  the  first  plot  in  this 
appendix  covers  this  special  case,  under  the  labeling  K  =  4, 

N  =  1,  M  -  1.  The  other  5  values  of  N,  along  with  the  3  possible 
combinations  of  K  and  M,  yield  15  additional  ROCs,  for  a  total  of 
16  ROCs  in  this  appendix. 
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860 


Figure  C-l.  FOCs  for  K  =  4,  N=l,  M  =  1  Figure  C-2.  ROCs  for  k 


Figure  C-5.  FOCs  for  K  =  4,  N  =  16,  M  =  1  Figure  C-6.  ROCs  for  K 
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Figure  C-9.  ROCs  for  K  =  2,  N=8,  M=2  Figure  C-10.  ROCs  for  K 
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Figure  C-13.  ROCs  for  K  =  l,  N  =  4,M  =  4  Figure  014.  ROCs  for  K 
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Figure  C-15.  ROCs  for  K  =  1,  N  =  16,  M  =  4  Figure  C-16.  ROCs  for  K 


APPENDIX  D  -  ROCS  FOR  KM  =  16,  PHASE-INCOHERENT  SIGNAL 


This  appendix  contains  the  ROCs  for  or-ing  with  pre-  and 
post-averaging  when  the  time-bandwidth  product  KM  is  fixed  at  16; 
the  possible  combinations  (from  table  1)  are  repeated  here: 


K 

M 

N 

16 

1 

1,2,4,8,16,32 

8 

2 

4 

4 

2 

8 

1 

16 

For  N  =  1,  only  the  product  KM  matters;  the  first  plot  in  this 
appendix  covers  this  special  case,  under  the  labeling  K  =  16, 

N  =  1,  M  «  1.  The  other  5  values  of  N,  along  with  the  5  possible 
combinations  of  K  and  M,  yield  25  additional  ROCs,  for  a  total  of 
26  ROCs  in  this  appendix. 
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Figure  D-l.  ROCs  for  K  =  16,  N  =  1,  M  =  1  Figure  D-2.  ROCs  for  K 
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Figure  D-5.  ROCs  for  K  =  16,  N  =  16,  M  =  1  Figure  D-6.  ROCs  for  K 
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Figure  D-17.  ROCs  for  K  =  2,  N  =  2,  M=8  Figure  D-18.  ROCs  for  K 
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Figure  D-21.  ROCs  for  K  =  2,  N  =  32,  M  =  8  Figure  D-22.  ROCs  for  K 
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Figure  D-25.  RDCs  for  K  =  1,  N  =  16,  M  =  16  Figure  D-26.  ROCs  for  K 


APPENDIX  E  -  ROCS  FOR  KM  *  64,  PHASE-INCOHERENT  SIGNAL 

This  appendix  contains  the  ROCs  for  or-ing  with  pre-  and 
post-averaging  when  the  time-bandwidth  product  KM  is  fixed  at  64; 
the  possible  combinations  (from  table  1)  are  repeated  here: 


K 

M 

N 

64 

1 

1,2,4,8,16,32 

32 

2 

16 

4 

8 

8 

4 

16 

2 

32 

1 

64 

For  N  =  1,  only  the  product  KM  matters;  the  first  plot  in  this 
appendix  covers  this  special  case,  under  the  labeling  K  *  64, 

N  =  1,  M  *  1.  The  other  5  values  of  N,  along  with  the  7  possible 
combinations  of  K  and  M,  yield  35  additional  ROCs,  for  a  total  of 
36  ROCs  in  this  appendix. 
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Figure  E-l.  ROCS  for  K  =  64,  N  =  1,  M  =  1  Figure  E-2.  ROCs  for  K 


Figure  E-3.  FOCs  for  K  =  64,  N  =  4,  M  =  1  Figure  E-4.  ROCs  for  K 


Figure  E-7.  ROCs  for  K  =  32,  N  =  2,  M  =  2  Figure  E-8.  ROCs  for  K 


Figure  E-9.  ROCs  for  K  =  32,  N  =  8,  M  =  2  Figure  E-10.  ROCs  for  K 


Figure  E-ll.  ROCs  for  K  =  32,  N  =  32,  M  =  2  Figure  E-12.  ROCs  for  K 


Figure  E-13.  FOCs  for  K  =  16,  N  =  4,  M  =  4  Figure  E-14.  ROCs  for  K 


Figure  E-15.  FOCs  for  K  =  16,  N  =  16,  M  =  4  Figure  E-16.  ROCs  for  K 


E-10 


Figure  E-17.  ROCs  for  K=8,  N  =  2,  M  =  8  Figure  E-18.  ROCs  for  K 
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Figure  E-21.  FOCs  for  K  =  8,  N  =  32,  M  =  8  Figure  E-22.  ROCs  for  K 
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Figure  E-25.  ROCs  for  K  =  4,  N  =  16,  M  =  16  Figure  E-26.  ROCs  for  K 


Figure  E-27.  ROCs  for  K  =  2,  N  =  2,  M  =  32  Figure  E-28.  ROCs  for  K 
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Figure  E-29.  ROCs  for  K  =  2,  N  =  8,  M  =  32  Figure  E-30.  ROCs  for  K 


Figure  E-33.  ROCs  for  K  =  1,  N  =  4,  M  =  64  Figure  E-34.  POCs  for  K 


E-19/ (E-20  blank) 


Figure  E-35.  HOCs  for  K  =  1,  N  =  16,  M  =  64  Figure  E-36.  ROCs  for  K 


APPENDIX  F  -  ROCS  FOR  KM  =  256,  PHASE- INCOHERENT  SIGNAL 


This  appendix  contains  the  ROCs  for  or-ing  with  pre-  and 
post-averaging  when  the  time-bandwidth  product  KM  is  fixed  at 
256;  the  possible  combinations  (from  table  1)  are  repeated  here: 


K 

M 

N 

256 

1 

1,2,4,8,16,32 

128 

2 

64 

4 

32 

8 

16 

16 

8 

32 

4 

64 

2 

128 

1 

256 

For  N  ■  1,  only  the  product  KM  matters;  the  first  plot  in  this 
appendix  covers  this  special  case,  under  the  labeling  K  =  256, 

N  «=  1,  M  -  1.  The  other  5  values  of  N,  along  with  the  9  possible 
combinations  of  K  and  M,  yield  45  additional  ROCs,  for  a  total  of 
46  ROCs  in  this  appendix. 
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Figure  F-3.  FOCs  for  K  =  256,  N  =  4,  M  ■  1  Figure  F-4.  ROCs  for  K  =  256 


Figure  F-5.  ROCs  for  K  =  256,  N  =  16,  M  =  1  Figure  F-6.  ROCs  for  K  =  256 
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Figure  F-7.  ROCs  for  K  =  128,  N  =  2,  M  -  2  Figure  F-8.  ROCs  for  K  =  128 
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Figure  F-9.  HOCs  f or  K  =  128,  N  =  8,  M  =  2  Figure  F-10.  KOCs  for  K  =  128 


Figure  F-13.  ROCs  for  K  =  64,  N  =  4,  M  =  4  Figure  F-14.  ROCs  for  K 
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Figure  F-17.  FOCs  for  K  =  32,  N  =  2,  M  =  8  Figure  F-18.  FOCs  for  K 


Figure  F-19.  ROCs  for  K  =  32,  N  =  8,  M  =  8  Figure  F-20.  ROCs  for  K 


Figure  F-23.  ROCs  for  K  =  16,  N  =  4,  M  =  16  Figure  F-24.  ROCs  for  K 
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Figure  F-25.  ROCs  for  K  =  16,  N  =  16#  M  =  16  Figure  F-26.  ROCs  for  K 
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Figure  F-29.  ROCs  far  K  =  8,  N  =  8,  M  =  32  Figure  F-30.  ROCs  far  K 
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Figure  F-35.  ROCs  for  K  =  4,  N  =  16,  M  =  64  Figure  F-36.  POCs  for  K 


660 


T-  lO 

p  o 
Q  P 


S  5 


to 

LLi 


F-37.  FOCs  for  K  =  2,  N  =  2,  M  =  128  Figure  F-38.  ROCs  for  K 


Figure  F-39.  ROCs  for  K  =  2,  N  =  8,  M  =  128  Figure  F-40.  POCs  for  K 


0.99 


F-22 


Figure  F-41.  ROCs  for  K  =  2,N  =  32,  Ms  128  Figure  F-42.  ROCs  for  K  =  1,  N  =  2,  M  =  256 
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Figure  F-45.  HOCs  for  K  =  1,  N  ■  16,  M  =  256  Figure  F-46.  ROCs  for  K 


APPENDIX  G  —  ACCURATE  EVALUATION  OF  A  CHARACTERISTIC  FUNCTION 
DIRECTLY  FROM  SAMPLES  OF  A  CUMULATIVE  DISTRIBUTION  FUNCTION 

Let  (real)  RV  x  have  (complex)  CF  f(£)  and  (real)  CDF 
c(u)  =  Prob(x  <  u).  The  only  quantities  available  are  samples 
c(nA)  for  all  n,  which  are  presumed  to  be  exact.  (The  PDF  p(u) 
of  x,  or  samples  of  the  PDF  p(u),  are  not  available.  Also, 
sample  values  (c(nA)}  are  not  taken  close  enough  to  yield  a 
decent  approximation  to  the  PDF  by  taking  differences.) 

Before  trying  to  relate  CF  f(S)  directly  to  the  available 
samples  of  CDF  c(u),  construct  a  model  CDF  cQ(u)  that  is  a 
reasonably  good  fit  to  c(u)  and  that  has  a  corresponding  closed 
form  CF  fQ(S);  in  addition,  CF  fQ(S)  should  decay  rapidly  to  zero 

as  l  ->  ±®.  Then,  the  desired  CF  at  real  argument  £  is  given  by 

f(£)  =  J  du  exp(i^u)  p(u)  =  J  du  exp(i^u)  c(u)  = 

=  J  du  exp(i^u)  ^j[cq(u)  +  c(u)  “  cQ(u)]  = 

«  fQ(S)  -  J  du  exp(iSu)  [c(u)  -  cq(u)]  ,  (G-l) 

using  integration  by  parts  and  the  fact  that  the  difference  of 
CDFs  tends  to  zero  as  u  ->  +®.  For  later  use,  define  the 
auxiliary  function 

f(S)  -  r 

g(S)  s  - -  =  -  J  du  exp(i£u)  [c(u)  -  cq(u)]  .  (G-2) 

At  the  origin,  g(£)  is  finite,  namely,  g(0)  =  /j  -  /jq,  which  is 

the  difference  of  the  means  of  the  respective  CFs  f(S)  and  f  (£). 
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SAMPLING  AND  ALIASING 


Define  an  approximation  to  the  desired  CF  f(5)  by  sampling 
the  integral  in  equation  (G-l)  at  increment  A  in  u  and  by  using 
the  trapezoidal  rule  (compare  equations  (B-14)  and  ( B— 1 7 ) ) : 

?(£)  s  fQ(5)  -  i£A£  exp(i£nA)  [c(nA)  -  cQ(nA)]  ,  (G-3) 

n 

where  the  sum  on  n  runs  over  +®.  (In  practice,  the  summation  on 
n  is  conducted  until  |c(nA)  -  cQ ( nA ) |  is  small  enough  that  the 
truncation  errors  on  both  tails  are  negligible.)  Expression 
(G-3)  can  be  manipulated  as  follows: 

2(5)  =  fQ(5)  -  i5  J  du  exp(i5u)  [c(u)  -  cQ(u)]  A  $A(u)  = 

=  f0(5)  +  i5  Y  g(5  -  n^fl  ,  (G-4) 

n 

where  $^(u)  is  an  infinite  train  of  delta  functions  at  spacing  A. 
This  equation  utilizes  the  fact  that  a  Fourier  transform  of  a 
product  is  the  convolution  of  the  corresponding  Fourier  trans¬ 
forms.  Also,  use  of  the  Fourier  transform  pair  (G-2)  was  made. 

When  the  explicit  representation  of  g(5)  in  equation  (G-2)  is 
used  for  the  n  =  0  term  in  equation  (G-4),  there  follows 

2(5)  =  fU)  +  U  £  gfe  -  n^fl  .  ( G-5 ) 

n*0  ^ 

All  the  manipulations  of  the  approximate  CF  2(5)*  from  its 
definition  (G-3)  through  equation  (G-5),  are  exact.  It  can  be 


immediately  seen  from  equation  (G-5)  that  the  summation  of 
aliasing  lobes  of  g(S)  for  n  /  0  constitutes  the  (complex)  error 
in  approximation  ?(£)  relative  to  its  desired  value  f(£).  To 
minimize  the  error,  function  g(£)  should  have  rapidly  decaying 
skirts,  so  that  the  aliasing  lobes  in  equation  (G-5)  do  not 
significantly  overlap  the  (-Jt/A,n/A)  region  of  £  where  f(£)  is 
substantial.  Recall  that  |f(0|  attains  its  maximum  value  of 
unity  at  the  origin. 

Function  fQ(£)  can  be  selected  as  desired.  As  a  starter,  it 
could  be  taken  as  a  Gaussian  CF,  with  its  mean  and  variance 
fairly  close  to  those  of  the  original  RV  x.  The  rates  of  decay 
of  |g(S)|  in  equation  (G-2)  and  of  |c(u)  -  cQ(u)  |  for  large 
arguments,  in  their  respective  domains,  would  then  probably  be 
dominated  by  the  specified  functions  in  this  problem,  namely, 

|f ( S) |/S  and  c(u) . 


COMPUTATION  OF  APPROXIMATE  CHARACTERISTIC 
FUNCTION  VIA  A  FAST  FOURIER  TRANSFORM 

Approximation  ?(£),  defined  by  equation  (G-3),  can  be 
calculated  at  whatever  £  values  are  desired.  However,  since  only 
the  interval  (-n/A,n/A)  in  £  is  of  interest  or  can  be  used, 
equation  (G-3)  is  evaluated  at  increment  A^  *  2n/(N1A)  according 
to 

=  fo(i^f)  "  £  exp(i2nmn/N1)  (c(nA)  -  cQ(nA)]  ,  (G-6) 
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where  integer  controls  the  spacing  of  the  samples  in  £.  Only 
the  values  for  |m|  <  N^/2  are  of  interest;  in  fact,  the  range 
0  <  m  <  N1/2  suffices  because  of  the  conjugate  symmetry  of  all 
the  functions  of  £  involved  in  equation  (G-3).  The  values  of  the 
summation  in  equation  (G-6),  for  this  m  range,  are  easily 
realized  by  means  of  an  N^-point  FFT  of  the  prealiased 
differences  of  CDFs  c(u)  and  c0(u);  this  prealiasing  feature  is 
incorporated  to  minimize  or  eliminate  truncation  errors.  Then, 
combination  of  the  FFT  values  with  the  sampled  fQ( O  function 
yields  the  approximation  £(*,)  to  the  desired  CF  f(S)  at  the 
values  £  =  mA^.  Since  the  FFT  in  equation  (G-6)  will  yield 
summation  values  for  0  <  m  <  N--1,  it  is  recommended  that  this 
entire  range  of  equation  (G-6)  be  plotted,  so  that  the  left  skirt 
of  the  first  aliasing  lobe,  namely,  n  =  1  of  equation  (G-5),  can 
be  observed  and  monitored. 


EXAMPLE 

An  example  with  a  Gaussian  CF  and  a  model  Gaussian  CF  with 
mismatched  parameter  values  was  used  as  a  test  case;  it  is 

fK)  =  exp(i /jS,  -  a2 12/ 2)  ,  c(u)  =  $  (U  ^ j  , 

fQ(S)  =  exp(i//QS  -  <t2E,2/2  )  ,  c0(xi)  =  — - — -)  .  (G-7) 
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The  plot  in  figure  G-l  reveals  the  desired  mainlobe  of  the 
magnitude  of  the  CF,  centered  at  the  origin,  with  magnitude 
errors  of  the  order  of  E-15  when  0  <  £  <  n/A.  For  larger  the 

first  aliasing  lobe  of  g(£)  dominates,  giving  useless  values  for 
the  approximation  f(£),  as  expected. 

To  support  the  claim  that  samples  {c(nA)J  of  the  CDF  may  not 
have  been  taken  fine  enough  to  achieve  a  decent  approximation  to 
the  PDF  p(u)  by  taking  differences,  the  CF  was  evaluated  (using 
the  trapezoidal  rule  for  integration)  for  the  PDF  estimate 
obtained  according  to  the  differences 

p((n-*s)A)  S  c(nA)  ~AC(  (n~1)A)  .  ( G-8 ) 

The  results,  displayed  in  figure  G-2,  indicate  a  very  significant 
error  between  the  approximate  CF  obtained  by  use  of  equation 
(G-8)  and  the  exact  CF  f(£).  The  maximum  magnitude  error  has 
increased  from  IE-15  to  3E-3,  which  is  more  than  10  orders  of 
magnitude.  The  utility  of  equations  (G-3)  and  (G-6)  is  made 
apparent  by  this  numerical  example. 

An  alternative  method  (to  that  in  this  appendix)  of  obtaining 
the  CF  directly  from  CDF  samples  is  afforded  by  the  modification 
in  equation  (B-20)  et  seq.  When  applied  to  example  (G-7)  here, 
excellent  results,  comparable  to  figure  G-l,  were  obtained. 
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APPENDIX  H  -  MATLAB  PROGRAM  FOR  EVALUATION  OF 


RECEIVER  OPERATING  CHARACTERISTICS 

The  numerical  procedure  utilized  is  based  upon  the  alias- 
suppression  technique  presented  in  appendix  B.  In  particular, 
smoothing  function  b^(u)  in  figure  B-l,  along  with  the 
corresponding  CF  relations  (B-14)  through  (B-17),  was  used  to 
evaluate  the  ROCs  in  appendixes  C  through  F. 

Every  figure  in  each  set  of  ROCs  has  three  parameter  values 
listed  on  it,  which  were  used  to  generate  those  specific  curves. 
They  are  the  bias  b,  sampling  increment  du,  and  FFT  size  n. 
Because  an  additive  constant  to  the  decision  variable  in  figure  1 
does  not  change  the  ROC  (reference  1),  bias  b  was  added  to  w(t) 
and  chosen  so  as  to  make  the  sum  variable  as  small  as  possible, 
but  never  negative  (within  roundoff  error).  This  procedure 
maximizes  the  unaliased  region  when  the  final  FFT  is  conducted 
from  the  CF  domain,  equation  (36),  to  the  EDF  domain  for  system 
output  w(t),  as  well  as  eases  the  requirements  on  du  and  n. 

Sampling  increment  du  is  applied  in  the  u  domain  to  equations 
(26)  and  (27);  it  must  be  taken  small  enough  that  the 
approximation  to  CF  fv($)  in  equation  (17)  is  not  aliased 
significantly  in  the  £  domain.  Finally,  FFT  size  n  must  be  taken 
large  enough  so  that  the  u-domain  span,  n  du,  in  the  final  FFT  of 
equation  (36)  from  the  CF  to  the  EDF  guarantees  insignificant 
aliasing.  Complete  details  on  this  cascade  procedure  are 
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presented  in  appendix  B  of  reference  1. 


All  three  parameter  values  have  been  chosen  by  a  trial  and 
error  procedure,  in  which  intermediate  CF  is  observed  in 

the  £  domain,  and  output  EDF  e  (u)  is  observed  in  the  u  domain. 

W 

Repeated  trials  led  to  the  tight  parameter  values  listed  on  each 
set  of  ROCs.  These  values  can  serve  as  starting  or  reference 
points  for  evaluation  of  additional  ROCs  with  different  values  of 
K,  N,  M,  and  p  that  are  of  interest  to  a  user.  The  following 
listing  of  the  MATLAB  program  corresponds  to  figure  E-36. 


H-2 


clear,  elf  %  MJWC  TR  11,166 

b=-190;  %  Additive  bias  to  w 

du=0 . 2 ;  %  Sampling  increment  in  u 

n=2^11;  %  EFT  size 

kc=l;  %  K,  amount  of  pre-averaging 

nc=32;  %  N,  amount  of  or-ing 

mc=64;  %  M,  amount  of  post-averaging 

rmin=-2;  %  Minimum.  SNR  (dB) 

rinc=.5;  %  Increment  in  SNR  (dB) 

nurm=16;  %  Number  of  ROCs 

disp(  'b  du  n;  kc  nc  me: ' ) 

disp(  [b  du  n;  kc  nc  me] ) 

xg=[le-6  le-5  le-4  .001  .002  .005  .01  .02... 

.05  .1  .2  .3  .4  .5]; 

yg=[le-6  le-5  le-4  .001  .002  .005  .01  .02... 

.05  .1  .2  .3  .4  .5  .6  .7  .8  .9  .95  .98  .99]; 
xg=phiircv(:xg)  ; 
yg=phiinv(yg)  ; 

[Xg,Yg]=meshgrid(xg,yg)  ; 

Magcf=zeros (n,  1) ; 

EDF=zeros  (n,  1) ; 

Bd=zeros  (n,nurrn-l)  ; 

nl=n-l;  n2=n/2;  kcl=kc-l;  ncl=nc-l;  n3=n2+l; 

pn=pi/n;  d2=du/2; 

arg=pn.* [I:n2] ' ; 

sinc=sin(arg)  ./arg; 

dispC  isnr  edfO  -  1') 

for  isnr=0:num 

rhodb=rrnin+rinc*  (isnr-1)  ; 
rho=10^  ( .  l*rhodb)  ,* 
if  (isnr=0)  rho=0;  end 
aa=sqrt  (2*kc*rho)  ; 

X=zeros  (n,  1)  +i*zeros  (n,  1)  ; 
meanv=0;  cdfk=0;  k=-l; 
while  (cdfk< .  5  |  area>le-20) 
k=k+l;  uk=du*k;  u2=uk+d2; 
sq=sqrt  (2*u2 ) ;  p0=es$)  ( -u2 ) ;  e0=p0 ; 
for  j=l:kcl 
p0=p0*u2/j ; 
e0=e0+p0; 
oxl 

c0=max(l-e0,0) ; 
if  (aa=0)  cl=c0; 
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else  [q  bes]=Qm(kc,aa,sq)  ; 

cl=max(l-q,  0)  ; 

axl 

cdfo=cdfk;  cdfk=cl*cOdicl; 
area=cdfk-cdf  o  ; 
j=mod(k,n)  ; 

X  ( j+1)  =X  ( j+1)  +area; 
meanv^earTR/'+area^Tjk;  %  uk,  not  u2 
end 

X=ff tgreen  (X)  ; 

Magcf=loglO  (X.  *conj  (X)  +le-50)  * .  5; 
plot  (Ffegcf) 

axis{[l  n+1  -16  0]);  grid  on 
paused) 

cbd=2*pi/  (n*du)  ; 

X(2:n3)=conj  (X(2:n3) )  ./sinc(l:n2) ; 
X(n3+l:n)  =0; 

X=X.^tfC; 

X(l:n3)=X(l:n3)  .*exp(i*b*dxi*  [0:n2]  '); 
meanw==meanv*rrc+b; 

X(1)=0; 

X(2:n3)=X(2:n3)  ./[l:n2] 

X=ff  tgreen  (X) ; 
a=.5-nneanw/  (n*du)  ; 
edf  0=a+litag  (X  ( 1 ) )  /pi  ; 
disp(  [isnr  edf 0-1] ) 
k=[l:n]'; 

edf=a-(k-l)  ./n+imag(X(k) )  ./pi; 

X=edf; 

EDF=loglO  (abs  (edf)  +le-30)  ; 
plot  (EDF) 

axis(  [1  n+1  -18  0] ) ;  grid  on 

pause  (1) 

edf=real(X(k) )  ; 

edf=mln  ( edf ,  l-le-12 )  ; 

edf =max  (edf ,  le-12 )  ; 

Bd( : ,  isnr+1)  =phiinv(edf )  ; 
axi 
be^p 
pause 
df 

bold  cn 

set(gcf,  'RaperPosition',  [.25  .25  8  10.5]) 


plot(xg,Yg,  'k' ) 
plot  (Xg,yg,  'k') 

plot(xg,xg,  'k')  %  zero  SNR  ROC 
plot(Bd(:  ,l),Bd(:,  [2:nurr+l] ) ,  'k') 
axis([xg(l)  xg(14)  yg(l)  yg<21)]) 
axis  off 
while  1 

T=input  ( 'pf  pd  isnr:  ' )  ; 
if  T(1)=0  break;  end 
pft=phiinv(T(l) ) ; 
pdt=phiinv  (T  (2 ) )  ; 
isnr=T(3)  ; 
for  j=0:nl 

if  (Pd(j+l,l)<pft)  break;  end 

93d 

xl=Pd(j+l/l)  ; 
x2=Pd(j,l) ; 
as=Bd  ( j  +1 ,  isnr+1 )  ; 
bs=Pd ( j +1 , isnr+2 ) ; 
cs=Ed  ( j ,  isnr+1 )  ; 
ss=  (cs-as)  /  (x2-xl)  ; 
al=pdt-ss*pf t ; 
fs=  (al+ss*xl-as)  /  (bs-as)  ; 
rc=rrrrin+rinc*  (isnr-l+fs)  ; 
disp([T  fs  rc]) 
and 

%  function  y=phiinv(x)  for  0  <  x  <  1 
%  y=l .  414213562373095*erf inv  (2*x-l)  ; 
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function  [q,bes]  =  Qn(m,a,b) 
start=le-100;  num=2000; 
a2=.5*a*a;  b2=.5*b*b; 
q2=start;  q3=start; 

for  j=l:ml 

q3=q3*b2/j; 

q2=q2+q3; 

arl 

ql=l;  q=q2;  bes=q3; 
for  j=l:num 
ql=ql*a2/j; 
q3=q3*b2/  ( j+ml)  ; 
q2=q2+q3 ; 
q4=ql*q2; 
q5=ql*q3; 
q=qK24; 
bes=bes4q5; 

if  (q4<=eps*q  &  q5<=eps*bes)  break;  end 
end 

q2=a2+b2 ; 
if(c£<708) 

el=exp  ( -q2 )  /start ; ' 
q=q*el; 
bes=bes*el; 
else 

ql=q2+log(start/q)  ; 

if  (ql>708)  q=0;  bes=0;  return;  end 

q=exp  (-ql)  ; 

ql=q2+log(  start /bes)  ; 

bes=e^>  (-ql) ; 
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